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Abstract. Let f! C K n be a bounded domain and for x £ O let t(x) be the 
expected exit time from f2 of a diffusing particle starting at x and advected 
by an incompressible flow u. We are interested in the question which flows 
maximize ||t||^oo(q), that is, they are most efficient in the creation of hotspots 
inside Q. Surprisingly, among all simply connected domains in two dimensions, 
the discs are the only ones for which the zero flow u = maximises ||t||^oo(qj. 
We also show that in any dimension, among all domains with a fixed volume 
and all incompressible flows on them, ||t||£,o° (n) ' s maximized by the zero flow 
on the ball. 



1. Introduction 

It is well-known that mixing by an incompressible flow enhances diffusion in 
many contexts. This is demonstrated, for instance, by the fact that the effective 
diffusivity of a periodic incompressible flow is always larger than diffusion in the 
absence of a flow [6], or that the principal eigenvalue fi u of the problem 



(1.1) 



-A6 + u ■ \76 = fi u 6 6 > in f£, 
6 = on dfl 



is never smaller than the corresponding eigenvalue no of (|1.1[) with u = 0. Classes 
of flows which are most effective in enhancing diffusion have been studied both on 
bounded and unbounded domains, and their characterizations have been provided 
in [HE]. 

On the other hand, it was observed in pQ that an incompressible flow may actu- 
ally slow down diffusion in the following sense. Consider the explosion problem 

-A6 + u-V6 = Xe't' infi, 
6 = on dfl. 

There exists A*(u) such that this problem has a solution for all A ^ A* (it) and no 
solution for A > A* (it) (see [5,9, 10 for u = and [5] for u ^ 0). Surprisingly, it was 
shown numerically in [T] that in a long rectangle there are incompressible flows with 
A* (it) < A*(0). This means that addition of a flow (which typically increases A* 
due to mixing) can sometimes instead promote the creation of hotspots and inhibit 
their interaction with the cold boundary dCl. 

The present paper is a step toward mathematical understanding of this diffusion 
slowdown effect of certain incompressible flows. We consider the problem 

-At" + u ■ Vr" = 1 in fl, 
t u = on dn 



(1.2) 



on a smooth bounded domain Q C R™, with u(x) an incompressible flow on U, (i.e., 
V • u = 0) which is tangential to dQ (i.e., u ■ n = on <9f2, with fi the outward 
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normal to 9fi). Physically, the solution t u {x) is the expected exit time from fi of 
the random process 

dX t = -u(X t ) dt + V2dB t , X = x, 

modeling the motion of a diffusing particle advected by the flow u. Although one 
might think that the expected exit time is always decreased by the addition of an 
incompressible flow due to improved mixing, this need not be the case. Our first 
result shows that in any bounded simply connected domain in K 2 which is not a 
disk, there are (regular) incompressible flows which increase the maximum of the 
expected exit time of X t from Q. 

Theorem 1.1. Let fl C K 2 be a bounded simply connected domain with a C 1 
boundary which is not a disk. Then there exists a C 1 , divergence free vector field 
v : f2 — » R 2 tangential to dfl such that Ht^IIl^q) > \\t ||j,<»(n)- 

Remark. We note that the incompressible flows are the natural class to study 
in this context. Indeed, if one considers general v (not necessarily divergence free), 
then it is easy to show that ||t"| can be made arbitrarily large by, for instance, 
taking v[x) — A(xq — x) with xq € Q and A sufficiently large. 

On a disk however, no (incompressible) stirring will increase this expected exit 
time beyond the one for u = 0. In fact we prove in any dimension that the L p -norm 
of the expected exit time can never be larger than that from a disk of equal volume 
with u = 0. 

Theorem 1.2. Let Q, C K™ be a bounded domain with a C 1 boundary and v : il — > 

R n a C 1 divergence free vector field tangential to dfl. Then for any p 6 [1, oo], 

|| «|| < \\ t o,d\\ 

\\ T \\LP(Q) ** ll T Wlp(D) 

where D C W l is a ball with the same Lebesgue measure as f2, and t° :D is the 
solution of (|1.3[) on D with u = 0. 

Remark. If D is a ball with Lebesgue measure V and center 0, then t°' D is given 
explicitly by the formula 



Ixl 2 



with r„ the Lebesgue measure of the unit ball in K n . 

There are, of course, other ways to quantify the effect of stirring on diffusion — 
see, for instance, [12] where many additional references can be found, especially to 
the physics literature. Closely related to the problem studied in the present paper is 
the following question. It is shown in [5] that for any p > n/2 there exists a constant 
C p (0) such that for any incompressible u tangential to d£l and any / € L p (fl), the 
solution of 

-A(p + u-Vcf> = f inn, 
= on dfl 



(1.3) 



satisfies ||0||l°° ^ C p (f2)||/||£p. It would be interesting to determine which flows 
achieve C p (f2) and how does C p (fi) depend on fi. Theorems 11.11 and 11.21 are a first 
step in this direction. 

The present paper is organized as follows. In Section [5] prove our main results, 
Theorems 11.11 and 11.21 Our proof of Theorem 11.11 involves a variational principle, 
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Proposition ^. 11 the proof of which is somewhat technical and therefore postponed 
to Section H This variational principle leads to an interesting PDE for the critical 
points of the expected exit time functional. We discuss properties of these critical 
points and provide some numerical examples in Section [3j 
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AZ was also supported by an Alfred P. Sloan Research Fellowship. GI was also 
partially supported by the Center for Nonlinear Analysis. 

2. Proofs of the main results 

2.1. Proof of Theorem 11.11 For a given incompressible C 3 flow u tangential to 
dQ, consider the family of Poisson problems 

-At Au + Au ■ Vt Au = 1 in Q, 

t Au = o on on 



(2.1) 



with A > 0. Let tp be the stream function of u, that is, i/> : f2 ^ K is C and such 
that ip(8Q) = and u = V^ip = {-d 2 ip,dnp). It is well-known OH] that if all 
critical points of ip are non-degenerate and no two of them lie on the same level set 
of ip, then the functions t Au converge uniformly to a limit f" which is constant on 
the level sets of ip and satisfies an asymptotic Freidlin problem on the Reeb graph 
of the function ip. If O is simply connected and ip has a single (non-degenerate) 
critical point (in which case either tp > or ip < on f2, and we will assume without 
loss the former), then we have the explicit formula 

r^iy) |o , , I 

(2.2) f"(y) = - r ^ - dh. 



( n Aipdx 

Here and elsewhere we let Sl^^h = {x <E £l \ tp(x) > h}, the /i-super-level set of 0. 
Notice that f u is just a reparametrization of if). As we prove in Proposition 13.11 
below, the formula (|2.2p holds also when the single critical point of tp is degenerate. 

We will start by considering only flows with the above property. That is, u is 
C 3 and such that the stream function ip only has a single critical point in fl (which 
is simply connected and ip > on Q). In particular, all super-level sets of ip are 
simply connected and ip attains a single maximum ip(xo) — M > 0. Moreover, any 
h E (0, M) is a regular value of ip and dil^p.h is a C 4 Jordan curve. 

Assume now that for some C 3 incompressible flow w the function t w has a single 
critical point and let ip = t w (which is C 4 ) and u = V^ip. Thus ip solves 

f -Aip + wVip = l infi, 
^ ' ' ip = on 90. 

Integrating this over il^ t h and using incompressibility of w, we obtain 

\^i>,h \ = Aipdx 



for any h £ (0, M). This together with (|2.2|) implies that f u = ip. That is, such 
solutions ip to the Poisson problem (|2.3p solve the Freidlin problem for themselves. 
We are particularly interested in the case w = 0, with ip = t° solving 

-At = 1 in A, 
t° = on dn 



(2.4) 
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and uq — V^t . Notice that then t° also solves 

-At + Au ■ Vr° = 1 in fl, 



(2.5) 



r° = on dfl. 



for any A £ ~R and so r° = f u ° . Let us therefore assume, for now, that is such 
that r° has a single critical point. 

We now assume that for any incompressible flow u on Q we have Ht^Hloc < 
||t°||l°o. In particular, ||f M ||_Loo < ||f n °||£oo for each u whose stream function has a 
single critical point. We will now show that this is the case only when O is a disc, 
thus proving Theorem II. II for all f2 such that t° has a single critical point. 

The key ingredient of our proof is that for all "infinite amplitude" expected exit 
times f", we have a variational principle which gives an explicit equation satisfied 
by the critical points (and thus the maximiser) of the functional I(ip) — ||f v ^||l°° 
(with ip having a single critical point). 

We set up the variational principle as follows. Let v (the "direction" of our 
variation) be any C 4 vector field tangential to dft. Let X be the flow (in the 
dynamical systems sense) given by 

^ = V oX e , X =Id. 
as 

Given a stream function ip with a single critical point, we perturb it by composing 
it with the flow X e . Let ip e = ipo X £ , u e = V^ip e , and f e = f u \ Notice that ip e is 
C and again has a single critical point (the maximum) Xq. Then f e also attains 
its maximum at Xq due to (|2.2p , so the variation of I in direction v is 

(2.6) V^,v)=^-f%xl) 

ae e=0 

We say that i/j is a critical point of / if for all C 4 (not necessarily divergence free) 
vector fields v tangential to <9f2, we have V(ip, v) = 0. Clearly any ip (with a single 
critical point) which maximises I is a critical point of /. So our aim is to prove 
that r° is not a critical point of / unless is a disc (assuming for now that r° has 
a single critical point). 

As mentioned earlier, the proof of this fact rests on obtaining an explicit equation 
for critical points of /. We can now do this by a direct computation using the 
Freidlin-Wentzel theory 0(5]. 

Proposition 2.1. Let fl C M 2 be a bounded simply connected domain with a C 1 
boundary and let ip > be a C 4 stream function on with tp(dil) — and a single 
critical point. Then ip is a critical point of the functional I if and only if 4> = f v ^ , 
the solution of the Freidlin problem (|2.2p with stream function ip, also solves 

-l 




(2.7) -2A#r) = l + |V0(x)| 2 / t^t, I / |V0| da 

We postpone the proof of Proposition 12.11 to Section [U but make two remarks 
before proceeding. 

Remark. One can also write down an explicit PDE (|4.4[) for the stream function ?p. 
This PDE, however, is somewhat more complicated, and we find it more convenient 
to work with (|2.7p involving the reparametrization </> of ip. 
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Remark. Assume that a C 3 flow w maximizes ||r' t "||ioo and r w has a single critical 
point. The argument following (|2.3p above then shows t w = f v r ™, so t w is a 
critical point of / and solves (|2.7p . 

Let now ip = r° have a single critical point xq € and assume that t/> is a 
critical point of /. Recall that ip = f u ° = f v ^ , so Proposition [2TT] implies that ^ 
solves (|2.7[) . Since — A-0 = 1, we obtain 

(2.8) |V^(x)| 2 / I^t f / IVVI derl = 1, 

immediately showing that |VV>| must be constant on the level sets of ip. Thus ^ 
solves the eikonal equation |V?/>(a;)| = g(ip) with g equal zero at the maximum of ip 
and positive elsewhere. It is well known that a solution of such equation does not 
have interior singularities only if is a disk and ip is radial [3] . In our situation this 
can be seen as follows. After reparametrization we may assume that g = 1, and ip 
attains its maximum at Xq . This introduces a singularity at Xq so let us suppose that 
ip does not have other interior singularities. Since the level sets of ip are connected, 
and the maximum is isolated, for any e > we can find a wavefront (a level set of 
ip) that is contained in a disc of radius e > centered at xo- By compactness, this 
wavefront is a positive distance e' € (0,e) away from xo- Absence of singularities 
now implies that we can evolve this level set, and the spheres of radius e and e' 
"outward" by the eikonal equation (with g = 1). Then each level set of ip obtained 
by this evolution lies entirely within distance s — s' < e from a circle. As e > is 
arbitrary we conclude that level sets of ip have to be circles. Since ip(dtt) — 0, we 
have that fl is a disk and ip radial. 

Thus we have proved that if t° has a single critical point, then it does not 
maximize / when O is not a disc. Since the claim of Theorem 1 1.1 1 for a disc follows 
from Theorem 11.21 (which we will prove shortly), we are left with considering the 
case of f2 such that r° has more than one critical point. We will use the following 
claim to reduce this to the previous case. 

Lemma 2.2. For any bounded simply connected f2 C K. 2 with a C 1 boundary, the 
set of maxima of t° is discrete. 

Proof. Let M — ||t°||l°c, let V — {x \ t°(x) = M}, and suppose V is not discrete. 
Since T> is positive distance from d£l, it has an accumulation point xo inside Q. 
Assume without loss of generality that a sequence x n € V converges to £o along 
the x-axis: (x n — Xo)/\x n — xo\ — > (1,0). Thus d x T°(x ) = d%T°(xo) = 0, so 
dyT°(xo) = 1 and the analytic implicit function theorem shows that there is a real 
analytic curve C containing x$ on which d y T° — (since r° is real analytic). It then 
follows that x n € C for all large n, and real analyticity of t°\c now shows C C T>. 

So T> contains an analytic curve which cannot end inside (by the previous 
argument) and must also stay away from dQ (by M > 0). This means that such 
a curve must be closed. But then r° > M inside the region enclosed by this curve 
(which is a subset of SI), a contradiction. □ 

We now return to the proof of Theorem 11.11 for general r° and assume that 
the zero flow maximizes ||t"||l°° among all incompressible flows u on Vt. We will 
reduce the problem to the previous case by showing that then the same is true for 
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a connected component of Qt-o^ containing a maximum of r , for all h sufficiently 
close to M = || r° || £00. We introduce some notation and make this precise below. 

Let t°(xq) = M for some x £ 12 and denote by 12 /j the connected component 
of £l T °,h containing xq. For any h and incompressible C 3 vector field w tangential 
to dQ. h , define Qn h ( w ) = umsu PA-»ooll T sirili 00 (nf,) (where t£™ satisfies (fL2|) with 
12 = Qh and u = Aw). Finally, choose ho < M sufficiently close to M so that 12/j 
contains no critical points of t° besides xq. 

Lemma 2.3. Assume that for all C 3 incompressible vector fields u tangential to 
912, we have ||t°||£oo(q) ^ Qn(u). Then for any C 3 incompressible vector field w 
tangential to dtt ha , we have ||r° - h \\L<»(n ho ) ^ Qa ho (w). 

Momentarily postponing the proof of Lemma l2~3l note that r° — ho is the expected 
exit time of Brownian motion, starting at x, from flh - That is, r° — ho is the 
solution of (|1.2|) with 12 = 12/j , and u — 0. Thus Theorem II. II for 12/j (which we 
have already proved) shows that 12^ is a disk of some radius, say R, and t° is 
radial in it. Since t°(x) + ^\x — Xo\ 2 is harmonic in 12 and radial near xq, it must 
be constant in 12. Since r°(912) = we have that 17 is a disk, completing the proof 
of Theorem 11.11 

It only remains to prove Proposition 12.11 and Lemma 12.31 Proposition 12.11 is 
proved in Section [¥J and we prove Lemma |2 . 31 below. 



Proof of Lemma \2.S[ The proof is based on the more general observation that 
changing any stream function near its maximum does not affect the asymptotic 
A — > 00 behavior of the solution of (|2.ip away from the maximum. We make this 
precise below. 

Let tfj be any C 4 function in 12 with -0(90) = and let xq, M, ho, 12/i be defined 
as above, with tp in place of r° (we will eventually choose ip = r ). For some 
hi £ (ho, M) let ip' be some C 4 function such that ip(x) = tp'(x) for x £ il\ 12/^, 
and denote u — V i/j, u' — V^tp' . Let ta and t' a solve 

, 2 9 x f -Ar A + Au ■ Vt a = 1 in 0, 

\ ta = on 912, 

and 
(2.10) 
We will first show 

(2.11) \\Vt a - Vr;|| L2(SAO;io) > as ^4 > +00, 

which, as mentioned earlier, says that perturbations of the stream function near xo 
do not affect the asymptotic A — > 00 behavior away from xq- 

To prove (|2.11|) . let <f>A = ta — t' a so that for any h £ [ho, hi] we have 



Au' ■ Vt' a = 1 in 12, 
4 = on 912. 



(2.12) 



—Acf)A + Au ■ S7(j>A = in 12 \ Q h , 
4>a = on 912, 

/ (V0 A -n)da = 0. 
Jdfi h 

where the third equation is obtained by integrating the difference of (|2.9p and 
(|2.10| over 12^, and using u ■ n = u! ■ n = on 912^. Multiplying (|2.12p by 0a and 
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integrating by parts, we obtain: 

f \V<j) A \ 2 = -f cj) A (V(j) A -n)da. 
Jn\n h Jon h 

Combining the flux condition in (|2.12[) with the last equality, we obtain: 



L 



I V(j)A | dx 



n\n h 



where 



1 



da 



dn h \dfl h \ Jda h 

is the streamline-averaged 4> A . Integrating this identity for h € [ho, hi] we obtain: 



(hi 



n\n h0 



a) 
in 

ho JdQh 



Qh \^h 1 



\\7<j) A \ z dx dh 



n\n h 



f>A — 4>a){ 1 ^4>a ■ n) da dh ^ 









PI 


&A - 4>A 


\V4>a\ da dh 


Iho Jan h 











2 \ 


I 


4>a - 4>A 


dx 


\Jn ho \n h 


i 





1/2 



Multiplying (|2.9p by ta, (|2.10p by t' a , and integrating over f2, we obtain the uniform 
bound ||Vryi||i,2(n), ||Vt^|| L 2 (0 ) C. Hence \\V (j) A \\ L ^(Q.) < C and it follows that 



(2.13) 



n\n ho 



^ C 



^a - 4>/ 



L^(n h0 \n hl ) 

< C (j|r A - TA|| L 2 (Oh0 \ nhi) + \\t' a - T A \\ L 2 (nho \ nhi) ^ 

We claim now that right side of (|2. 13(1 tends to zero as A — > oo. Indeed, multiplying 
(|2.9|) by u ■ Vta, integrating, using incompressibility of u, and the fact that u ■ n = 
on c?f7 gives 

A [ (u-Vr A ) 2 dx= f (u ■ S/T A )Ar A dx — — f -^—^ — (u-Vr A )dx 
Jn Jn Jn oxj oxi 



8t a du m dr A 
n dxj dxj dx m 



dx^C [ \Vt a \ 2 sC C. 
Jn 



As \u\ is strictly positive in Qh \ &hn it follows that 
( 2 - 14 ) \\ T A-r A \\ L 2 {nho \ nhi) -> 

as A — ► +oo. The argument for t' a is identical, completing the proof of (|2.11D . 

In order to improve the iJ 1 (f2 \ Qh ) bound (|2. 1 1[) to a bound in L°°(Q, \ fi/ l0 ) 
we simply note that, given any e > and A > A n , using (|2.14|) we may find a 
streamline dflh' with h' < h but arbitrarily close to h^ so that 

\\t A - T A \\mdn h ,) + \\t' a ~ T A \\ L 2 {nhl) + ||Vta - Vr^||i2(gn h ,) < e. 

It follows that then 



\t a - T A \\L<*>(n h ,) + \\t' a - ^AlU~(n h ,) < Ce, 
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and, in addition, \ta — f A \gn h , < Ce because of (|2.11j) and since ta = t' a = 
on dfl. Finally, since ta and t' a satisfy the same equation outside of Q \ Qh , the 
maximum principle implies that \ta — t' a \ < Ce in f2 \ flh - 

Now assume that tp — t° maximizes ||t u ||£oc (then uq = V^t maximizes Qn) 
but wo is not a critical point of /. Then there exists a C 4 stream function ip' on f2, 
equal to t° on £1 \ fi/j , such that for w = V^ip' (when restricted to £lh ), 

(2-15) M = \\r°\\ Lao <h +Qn ho (w). 

We can assume that t/j' has a single critical point in tth because so does r° as well 
as all the perturbations ip £ considered in the proof of Proposition ^. 11 Moreover, we 
can assume ip' = t° on fl \ for some hi > ho because it is sufficient to consider 
such perturbations in that proof (see the remark after the proof of Lemma . 

Then the previous argument shows \\t Aw — t° \ \ (n hQ _ s \n ho+s ) — * as A — > +oo. 
In particular, t Aw > h — S on d£lh f° r an large A, with 5 = (ho + Qn h (w) — M)/2. 
This means that t Aw > ho — 5 + Tq w on Qh by the maximum principle. But then 

M > Q n (w) >h -S + Qn h0 (w) = M + S> M, 
a contradiction. This finishes the proof. □ 

2.2. Proof of Theorem 11.21 We can assume that v is sufficiently smooth (and 
approximate general v with smooth ones). Let us denote r = t v and flh = ^r.h- 
Then by Sard's theorem the set A of regular values of r has full measure. Thus 
dflh is a finite union of sufficiently smooth compact manifolds without boundary 
for each h S A (moreover, A is then open because r S C 2 (f2)). 

Let fi* and r* be the symmetric rearrangements of f2 and r. That is, fi* is 
the ball with volume |f2| = V centered at the origin and t* : f2* — > i? + is the 
non-increasing radial function such that the ball S!£ = {a; € £1 1 t*(x) > /i} satisfies 
|fi? | = |f^| for each /i€M (with fi^ as above). 

Let now h € A. The isoperimetric inequality gives 

(2.16) \dn* h \^\dn h \, 

with equality precisely when is a ball. Since u is divergence-free and r is constant 
on dflh, we have 

(2.17) / |Vr|<f<7 = - / ^-da= j (-At + v-S/r)dx = \n h \ = \n* h \ . 



dQ. h Jon 



v 



Finally, the co-area formula yields 



_ 9Slh ,Vr| dh> 8h> J dn , |W. 

Thus by (|2.16[) and the Schwarz inequality, 

/ |Vr*| da [ -±-da = \dn* h \ 2 < |9fi fe | 2 < / |Vr| da [ da. 

In view of (|2.17p and (|2.18j) we obtain 

/ |Vr*| da^ [ \Vt\ da=\n* h \, 

JdQ.* h JdUh 

with equality precisely when Q % is a ball and |Vr| = —dr/dv is constant on <9f2/j. 
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So if j(\x\) = t*(x) and p = (V/T n f/ n is the radius of ft*, then with £„ = nT n 
the surface of the unit sphere, 

F r n r 

7 (p) - and < - 7 '(r) < = - 

when 7(r) G .4. Since A has full measure and 7 is continuous, we have 

2 2 

(2.19) 7 (r) ^ P -^- = 7(r) 

for all r € [0, p], with 7 = 7 precisely when all Sl^ are balls and r is radial (thus so 
is v(x) ■ x, hence v(x) ■ x = since v is divergence- free) . Now (I2.19[) gives 

|n h | = \n* h \ = \{x e n* 1 7 (M) > Ml < K* e ^* 1 7(M) > Ml 

and the claim follows. □ 

3. Properties of the Maximizer 

We start by proving (|2.2p . 

Proposition 3.1. Let ip > be a C 4 stream function on a bounded simply con- 
nected domain f2 C K 2 luzi/i a single critical point and let u = V^V- Then t Au — > f u 
uniformly on Q, where f w is given by 

(3 - i) ™--L j5^ dh - 

Proof. Assume first that the maximum of 1^ is non-degenerate and let sup ip = M > 
0. It is then proved in [2] that, as A — + 00, the functions t Au converge uniformly 
on £1 to f u with f u (y) — ?(ip(y)), where f solves the effective problem 

1 "US) -1, 



(3.2) { T(h) dh V dh / 
f (0) = and f is bounded on (0, M) 

on the interval (0, M), with the coefheients 

(3.3) T(M= / p^T, p(fc) = / \Vi>\da. 
By Green's formula and 



p(h) = \Vip\da = - Vip-hda = - Aipdx, 

where we have used V^> • n = — \Vip\ on dCl^^. By the co-area formula, 

T(r) dr = dx, 



and thus (|3.2p reduces to 

(3.4) f > {h) = l™ T ^ dr + C - +C 



p{h) -f^A^dx' 

Non-degeneracy of the maximum of tp shows that h 1 M f Q t Atp dx stays bounded 
away from zero and infinity as h f M. Boundedness of f then forces C = 0, 
completing the proof of the non-degenerate case. 
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If the maximum of ip is degenerate, we let ipn be a C 4 stream function with a 
single non-degenerate critical point which agrees with ip on f2 \ fi,/,.M-(i/n)- The 
proof of Lemma l2~3l with u = S7 ± ip and v! = u n = V" L, I i, shows that T Au — T Au ^ _^ 
as A —>■ oo, uniformly on Cl$ t M-(2/n)- But 



/ Atpdx = - / \Vtp\da 



shows that f u and r" n coincide on Sl$ t M-(i/n)> so rA " —> t v as A — > oo uniformly 
on Slip t M-(2/n)- The result now follows by taking n — > oo and noticing that for 
large n and large A, the oscillation of t Au on n^,M-(2/n) nas to be small thanks 
to the small oscillation of t Au on 9f2^,M-(2/n)) small diameter of fi^ Af_(2/ n )) and 
the maximum principle. □ 



We are presently unable to analytically prove existence of solutions to (|2.7|) . The 
structure of the nonlinear term in (|2.7p yields itself naturally to some apriori esti- 
mates. These, however, are not strong enough to prove existence, mainly because 
they do not seem to provide any form of compactness. 



Proposition 3.2. Let <fi be a C solution of (|2.7| with a single critical point and 
= o?i dfl. Then 

(1) Hh^§. 

(2) For any Borel function f , 

f{4>) \Acj>\ dx= [ f(<f>)dx, 



and, in particular, I |A0| dx = \Cl\. 
Jn 

(3) If t satisfies — At = 1 in Q with r = on dfl, then 

/ mega \A(p - Ar| < |0| . 
Jo 

These estimates do give us some insight as to the nature of classical solutions 
to (|2.7p . For instance, the first two assertions give L°°(SY) and bounds on 

<j>, while the third is an explicit upper bound on the distance between a classical 
solution of (|2.7[) and the exit time of the Brownian motion from 57. 

Proof. The second assertion follows by multiplying (]2 . T[) by f{<f) and using the 
co-area formula. As a consequence, for any h > we have the identity 

\n^ h \ = / \Vc/)\ da. 

Then (1) follows by a rearrangement argument as in the proof of Theorem ll.21 For 
the last claim, note that 

2(At(») - A^)) = ' W0E)|2 / j^da-1. 



da 



|V0| 



By the co-area formula, the integral over fl of the first term is exactly Since 
that term is non-negative, the strict inequality in (3) follows. □ 
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Since an analytical proof of existence of solutions for (|2.7|) is at present intangi- 
ble, we turn our attention to numerics. As boundary integrals are problematic to 
compute numerically it is more convenient to work with the equation 

(3.5) - 2A<f>(x) = 1 + V<Kar) ■ V In | % (x) \ , 

which is equivalent to (|2.7p . Surprisingly an iteration scheme of the form 



-A0o = 1, 



= 



-2Acf) n+1 (x) = 1 + V<p n (x) • Vln|{0„ > <t>n(x)}\ , 4> 



n+l 



on 







does not always converge. For certain domains, it turns out that numerically 
ll^nlli 00 -*ooasn-» oo, which is clearly not representative of the solution of (|2.7|) 
as it violates the last assertion in Proposition 13.21 
It turns out that an iteration scheme of the form 



(3.6) 
(3.7) 

(3.8) 



-Ac 



1. 







2Acj> n+ x{x) = 1 + Vcjy n {x) • Vln|{0„ ^ <t>n(x)}\ , <A 



n+i 



"n+l 



{^„+i ^ h} 



i_— ny on 



dh. 



does converge rapidly to a numerical solution of ()3.5j) . In fact, 

by 



3.8j) can be replaced 



-Ac 



'n+l 



V0 



•n+l 



1 



for some large, fixed A, which produces better numerical results. Figure [3] shows 
contour plots of the solution to (|2.7p in two different domains. For comparison, the 
expected exit time from the domain tq is shown alongside each plot of <f>. 

We are unable to prove convergence of these numerical schemes, just as we can 
not establish existence of solutions of (|2.7p . However, one immediate observation 
from Figure[3]is that the level sets of <f> become circular near the maximum. Indeed, 
for any classical solution of (|2.7p . this must be the case. 

Proposition 3.3. Let <p be a smooth solution of (|2.7p and assume that (f> attains 
a local maximum at (0,0), then d xx (j)(0) = d yy 4>(0). 

Proof. We will show if 4> is any smooth function which attains a maximum at 0, 
then the last term in ()3.5[) is continuous near if and only if d xx cf>(0) = d yy (j)(Q). 
This immediately implies the proposition. 

Assume first that the Hessian of <fi at is not degenerate (in this case assuming 
4> £ C 3 will be enough for the proof). We rotate our coordinate frame and assume 
without loss of generality that 



4>{x,y) =M 



x 2 V 2 , , 
--- +C3 (x,y) 



where 03(2;, y) is some function involving only third order or higher terms. Now for 
any e > 0, define f(e) by 



(3.9) 



f(e) = \n*,M- e \ = 



+ ^2 ^£ + c 3 (^,y) 



irabe + 0(e 2 ), 
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(c) Maximiser ip (d) Expected exit time tq 



Figure 1. Maximisers and expected exit times from two different domains. 



whence 

(3-10) ±\n{f{e)) = \ + 0{l). 

Thus 



v<Kx, v ). vin|n^ (Biy) | = \v<Kx, y )\ 2 ( j,^ +Q(i; 



|V ^ y)|2 ■0(l)|V*(»,„)| a 



M-(j>{x,y) 

The second term on the right is certainly continuous at 0, since V</>(0) = 0. The 
first term is continuous at if and only if a = b. 

In the case that the Hessian of <f> is degenerate at 0, the above proof works 
with minor modifications. Using a higher order Taylor approximation of </>, the 
right hand side of (|3.9|) becomes c\e C2 for two constants c\ > 0, and C2 = \ + ^ 
with n 3. Now replacing 1/e with 02/2 in (|3.10[) , the remainder of the proof is 
unchanged. □ 

4. Proof of Proposition 12.11 
First, we obtain an expression for V(ip,v). Let Qh = f2^> ft and = fl^ 



Lemma 4.1. Let M = supip, then the variation (|2.6p is 
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(4.1) Vty,v) = 



A I 



° \fdn h at d a 



\n h \ J (Jr (v ■ V-0) - Aipv • hj da+ 



dip 

Tpr da 

an h on 



v ■ hda 



dn h 



dh. 



Proof. This follows from Proposition 13. II Note that 



< 4 - 2 > 5 



e=0 



m\ = — 

de 







[ dx = -r 


I detVX" 1 


e= o Jn h 



V • vdx = — J v ■ n da, 
ldn h 



and 



d 
de 



=o Jen% 



dip 5 
dh 



da = — 



=o -in 



Aip e dx 

/ (Ai/> £ ) o X~ l |det VX~ X | cfa; 

=0 "'fife 

A (w • W) - v ■ VAip - (V • v) Aip 

[ — (v-Vip)da — ( (Aip)v-hda 
Jan h on J drih 



dx 



(4.3) 



Thus, using (|3.1[) . and equations (|4.2p — (|4.3[) we are done. □ 

Before proving Proposition ^. 11 we require a lemma. 

Lemma 4.2. A C A stream function tp (with a single critical point) is a critical 
point of the functional I if and only if it solves 

(4.4) VF^ • + 2i^A</> - G 4 , = 0, 

where F^p and G^p are defined by 

Gip(x) = | / an ^ F i>( x ) = \^4>{x)\G4>(x) 2 - 

Proof. With Fy,G^, as above, equation (|4. 1 [) reduces to 
(4.5) 

" d 

— (v ■ V?/>) — Aip v ■ h 



V{tl>,v) = / F$ 



M 



dn h 



dh 







da dh + 


!" G " S 



By the co-area formula, we have, first, 



,1/ 



(4.6) J G^ J vhdadh 
second, 



M 



dn h 



G^, v-Vip dx, 



(4.7) 



M 



F^ J [-Aipv ■ h] dadh= \ F^ Aip (v ■ Vip) dx, 

dfl h 
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and, finally, 



(4.8) 



j J -^r(v ■ V^) dadh = - J^F^y(v ■ V^) • Vip dx 

(v ■ V-0) (V-F^ • VV^ + Fj,Aip) dx. 



In the last equality we used the identity 

(v ■ Vip) {VF^p • Vt/> + F^Aip) + F^p V(v • V^) • S/ip = V ■ [F$ (v ■ W) V^] 

and the fact that F^, (v ■ Vip) Vip — on 9fi. 
Using (|4.6p - (|4.8p . expression (|4.5[) becomes 

V(^», u) = / (« • W>) [VF V , • Vif) + 2F^,Aip - G$] dx. 



Thus V(ip,v) = for all C 4 functions i; which vanish on dil if and only if equa- 
tion (TO]) holds. □ 

Notice that the same conclusion is obtained if we ask V(ip,v) — only for all v 
compactly supported inside O. 

Proof of Provosition [2J\ Note that the variation V(ip, v) in (!2-6[) depends only on 
the geometry of the level sets of the stream function ip and thus is invariant under 
reparametrizations. Thus if ip is a solution of (|4.4[) , then for any monotone function 
/, foip is also a solution of (|4.4|) . Note that = f v ^, the solution of the Freidlin 
problem (|2-2[) with stream function ip, is only a reparametrization of the level sets 
of ip. Thus to prove Proposition 12.11 we only need to show that if ip solves (|4.4I) . 
then <f> solves (|2.7[) . 

Note that since cf> solves the Freidlin problem (|2.2[) . we have 



(4.9) 



— 



and so Fa, = — Ga,. We also have 



VF (z) = -VG^z) = G (a;) 2 V 

= -G^x) 2 v |n^(a;)| = 

= -G (x) 2 V0(a;) 
and using this in (I4.4[) immediately yields (|2.7|) . 




■ da dh 



□ 



We remark that any solution to (|2.7p is automatically a solution to the Freidlin 
problem with itself as stream function (i.e. satisfies (|4.9[) ). Indeed, integrating (|2.7p 
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over £l<t>,h and using the co-area formula gives 



M 

2 



-l 

da f da 



-2 j A =|««.i + y ki j w m j m \J m*\ *. 

and hence 



2 y iv0i=|o Mo |+y y -2-^^=21^1, 

showing (|4.9|) is satisfied. 
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